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Abstract 

The bright matter wave soliton propagation through a barrier with a rapidly oscillating position 
is investigated. The averaged over rapid oscillations Gross-Pitaevskii (GP) equation is derived. It 
is shown that the soliton is dynamically trapped by the effective double-barrier. The analytical 
predictions for the soliton effective dynamics is confirmed by the numerical simulations of the full 
GP equation. 
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I. INTRODUCTION 



The tunneling of quantum particles is one of the fundamental problems of quantum 
mechanics In connection with this problem, the propagation of matter wave packets 
through different types of potentials has recently attracted a lot of attention. In partic- 
ular it is interesting to study the transmission through time-dependent barriers. Indeed, 
experiments with dynamical tunneling of cold atoms in time-dependent barriers of the mul- 
tiplicative form f(t)V(x) show the possibility for the tunneling period control 

as. 



theory is developed for linear wavepackets in Ref. 



a. 



The 

Soliton propagation through such 
a barrier is investigated in Ref. jsj]. Another interesting problem is to consider a barrier 
with a rapidly oscillating position, i.e. a potential of the form V(x — f(t)). In the linear 
regime of propagation resonant effects are exhibited in the transmission process through an 
oscillating barrier It is therefore relevant to investigate the propagation of a nonlinear 

wave packet through an oscillating barrier. We consider this problem for a bright matter 
wave soliton in an attractive Bose-Einstein condensate (BEC) in the presence of a barrier 
potential with an oscillating position. Such a barrier can be achieved by using a laser beam 
with a blue-detuned far-off-resonant frequency. This generates a repulsive potential and the 
laser sheet plays the role of a mirror for cold atoms. Such a laser sheet has been used to 
study the bouncing of a BEC cloud off a mirror |8|. By moving the position of the mirror 
we obtain an oscillating barrier. 

In this work we will study the dynamical trapping of soliton in BEC with a rapidly 
oscillating barrier potential. The averaged over rapid oscillations Gross-Pitaevskii (GP) 
equation will be derived. Applying a variational approach to this equation, we will analyze 
the conditions for soliton trapping by the effective potential. The analytical predictions are 
checked by numerical simulations of the full GP equation. 

II. THE AVERAGED GROSS-PITAEVSKII EQUATION 

The BEC wavefunction ip in a quasi one-dimensional geometry is described by the GP 
equation 

m T + ^xx - V P (X - f p (T), T)V - gi D \^ = 0, 
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where V P (X,T) is a barrier potential. Here g\o = 2ha s u± is the one-dimensional mean 
field nonlinearity constant, a s is the atomic scattering length, uj± is the transverse oscillator 
period, and a± = \/H/(muj±) is the transverse oscillator length. Below we consider the 
case of BEC with attractive interactions, i.e. a s < 0. To avoid collapse in the attractive 
BEC, the condition \a s \N/a± < 0.67 should be satisfied, with iV the number of atoms. 
The barrier potential has an oscillating position described by the zero-mean time-periodic 
function f p (T) with period T p . We consider the case of a fast-moving potential in the sense 
that e := uj±T p -C 1. We also allow the barrier potential profile V P (X,T) to vary slowly in 
time, i.e. with a rate smaller than u±. Introducing the standard dimensionless variables 

V2X , m f —. i ^ T y v 



x = , t = uj±T, u = y/\a s \ip, V 



we obtain the GP equation with fast moving potential 

t 



iu t + u xx - V \ x - f(-),tj u + 2\u\ u = 0, (1) 

where _ 

fit) = —f P (T p t) 

CL± 

is a periodic function with period 1 and mean 0. The small parameter e describes the fast 
oscillation period in the dimensionless variables. We look for the solution in the form 

u{x, t) = w (0) O, t, -) + eu {1) (x, t, -) + £ 2 u {2 \x, t, -) + • • • (2) 

where u^°\ u^ l \ . . . are periodic in the argument r = t/e. We substitute this ansatz into 
Eq. ((J) and collect the terms with the same powers of e. We obtain the hierarchy of equations: 

iu® = 0, (3) 
iuf ] + uf x -V(x- f(r),t) u<® + 2| M (°)|V°) = -iu?\ (4) 
iu? + -V(x- f(r),t) + A{u^) 2 W) + 2\u^\ 2 u^ = -iu {2) . (5) 

The first equation (j3J) shows that the leading-order term does not depend on r. The 
second equation (j3J) gives the compatibility equation for the existence of the expansion (J2J): 

iuf ] + 4°2 - (V(x - f(-),t)) + 2\u^\ 2 u^ = 0, 

where (•) stands for an average in r. Denoting 

V cS (x)= [ V(x-f(r),t)dr (6) 
Jo 

3 



we obtain the averaged GP equation 

iu t + u xx — V c fi(x, t)u + 2|u| 2 m = 0. (7) 

The equation (jSJ) allows us to compute the first order correction, whose amplitude is of order 
e. 

If, for example, /(r) = asin(27rr), then a straightforward calculation shows that 

Vtf(x,t)=K eS *V(x,t), (8) 

where * stands for the convolution product (in x) and K e s is given by 

^eff(x) = \ 7 =L= l ( _ a>a) (a;). (9) 
7T Va 2 - x 2 

If, moreover, V(a;, t) is a delta-like potential centered at a position xo(t) that moves slowly 
in time, then 

V eS (x,t) = K eS (x - x (t)) 

is a moving double-barrier potential. By this way, one can generate a moving trap potential 
to manage the soliton position. 

III. VARIATIONAL APPROACH FOR THE SOLITON MOTION 

In absence of effective potential V e s = the NLS equation (J2J) supports soliton solutions 
of the form 

e i2fi(x-x s (t))+4i(^ 2 +u 2 )t 

Uo{x,t) = 2U cosh[2z/(x - x 8 (t))} 
with x s (t) = x s (0) + A/xt. In the presence of a stationary potential V e s(x), the soliton 

dynamics can be studied by perturbation techniques. Applying the first-order perturbed 

Inverse Scattering Transform theory [llj], we obtain the system of equations for the soliton 

amplitude and velocity 

dv du 1 . . , dx s 



where the prime stands for a derivative with respect to x and W v has the form 

f°° 1 z 
W v {x) = v I — -j— V eS (x - — )dz 
J-oo cosh (z) Iv 



= K u *V eS (x), (10) 
cosh {Zisx) 
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Therefore, we can write the effective equation describing the dynamics of the soliton center 
as a quasi-particle moving in the effective potential W UQ : 

^ = -<(*.)> (I 2 ) 

where 4i/ is the mass (number of atoms) of the initial soliton. This system has the integral 
of motion 

y {^j + W U0 {x s ) = ^l (13) 

where 4/j, is the velocity of the initial soliton. Note that this approach gives the same result 
as the adiabatic perturbation theory for solitons that is a first-order method as well. This 
adiabatic perturbation theory was originally introduced for optical solitons Q and it 
was recently applied to matter wave solitons 

The quasi-particle effective potential W u is given by 

W v (x) = (K u *K eS )*V(x). 

In the case where /(r) = a sin(27rr), the kernel K e g is given by and the Fourier transform 

of the quasi-particle potential has the explicit form 

TrkJo(ka) 
2sinh(g) 

where V is the Fourier transform of the reference potential V. The effective quasi-particle 
potential is plotted in Fig. [T] for a Gaussian barrier potential V. The potential W v has a 
local minimum at x = between two global maxima that are close to x — ±a. The trap 
amplitude is 

AW V = maxW„(x) - W u (0). 

X 

Let us consider an input soliton at x = with parameters (u ,fi ). If the initial soliton 
velocity is large enough, then the soliton escapes the trap. There is a critical value \x c for 
the initial soliton velocity parameter /i defined by 

that determines the type of motion: 

- If \fio\ < fi c then the soliton is trapped. Its motion is oscillatory between the positions ±Xf 
defined by 

W n (x,)-W n (0) = &$i4 (14) 





FIG. 1: Effective potentials V e s{x) for the averaged equation (J7| and W u (x) for the quasi-particle 
motion <(TU). Here v = 0.25, V(x) = 10 exp(-x 2 /4), and f(r) = 10sin(2vrr). 

- If |/!o| > /i c , then the soliton motion is unbounded. It escapes the trap and reaches the 
asymptotic velocity parameter /i a given by 

^o(O) 



iis 



(15) 



which shows that the transmitted soliton velocity is larger than the initial soliton velocity. 

As we shall see in the numerical simulations, if the initial soliton parameters are close to 
the critical case |/io| ~ /i c , then radiation effects become non-negligible. The construction of 
an efficient trap requires to generate a barrier potential V that is high enough so that AW„ 
is significantly larger than 8z/q/i . 

Finally, if the potential V is not stationary but has an explicit time- dependence, such 
as a drift, then the adiabatic equations derived in this section still hold true if the time- 
dependence is slow enough. 



IV. NUMERICAL SIMULATIONS AND EXPERIMENTAL PREDICTIONS 

In this section we show the motion of a soliton with initial parameters uq = /io = 0.25 (the 
soliton velocity is 1). The reference barrier potential is V(x) = 10 exp(— x 2 /4). Note that 
in this case, AW^ ~ 0.323 and \i c ~ 0.4. Therefore, the quasi-particle approach predicts 
trapping. 

In Fig. |21 the modulated potential is V[x — 10sin(50t)] (stationary trap). Here e = 
27r/50 ~ 0.126 and we check by full numerical simulations of the GP equation (JIJ that e is 
small enough to ensure the validity of the theoretical predictions based on the averaged GP 
equation (JJJ). We plot the soliton profile and the soliton mass, which shows that radiative 
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FIG. 2: Soliton profile in a narrow trap V[x — 10sin(50i)] (left picture). Mass decay (right picture). 
Bottom picture: Soliton motion predicted by the quasi-particle approach (dashed) and observed 
in the numerical simulations (solid). 

effects are very small. We also compare the soliton motion obtained by the numerical 
simulations with the one predicted by the quasi-particle approach, which shows very good 
agreement. 

In Fig. |3]the modulated potential is V[x — 10 sin(50t) — 0.25i] (moving trap, with velocity 
0.25). We plot the soliton profile and the soliton mass, which shows that radiative effects 
are very small. 

In Fig.|3]the modulated potential is V[x — 10 sin(50i) + 0.25£] (moving trap, with velocity 
—0.25)). Radiation effects are small, but not completely negligible. It can be seen that 
some radiation is transmitted through the edges of the effective double-barrier trap. We 
then double the potential amplitude and consider 2V[x — 10sin(50t) + 0.25t]. The new 
simulation is shown in Fig. where radiation effects are seen to be negligible. 

In Fig. IHlwe show the motion of a soliton with initial parameters v$ = 0.25 and fio — 0.5 
(the soliton velocity is 2). Since we have AW Uo ~ 0.323 and /i c = 0.4, the quasi-particle 
approach predicts unbounded motion. The results of the numerical simulations confirm this 
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FIG. 3: Soliton profile in a narrow moving trap V[x — 10sin(50t) — 0.25t] (left picture). Mass 
decay (right picture). 
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FIG. 4: Soliton profile in the moving trap V[x — 10sin(50i) + 0.25i] (left picture). Mass decay 
(right picture). 

prediction. We can observe that half the initial soliton mass is transmitted through the 
effective barrier during the first interaction, and the transmitted soliton has a velocity that 
is larger than 4/j, . 

Let us estimate the typical values of the parameters for an experimental configuration. 
We consider the case of a condensate with 7 Li atoms in the a cigar-shaped trap potential 
with uj± = 2-k x 1 kHz. In this case, e = 2n/50 corresponds to 1/T P = 50 kHz. The barrier 
height is V = lOhw^. The number of atoms can be taken as iV ~ 10 3 for the scattering 
length a s = —0.2 nm. The transverse oscillator length is a± — 1.2 /im, so the barrier width 
is d ~ 2 /im and the oscillations amplitude of the barrier position is L ~ 10 /xm. The velocity 
is normalized by the sound velocity which is c ~ 4 mm/s, so the soliton velocity parameter 
fj,o = 0.25 corresponds to a velocity equal to c. 
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FIG. 5: Soliton profile in the moving trap 2V[x — 10sin(50t) + 0.25t] (left picture). Mass decay 
(right picture). 
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FIG. 6: Soliton profile in the stationary trap V[x — 10sin(50i)] (left picture). Mass decay (right 
picture). Here the initial velocity parameter is /Iq = 0.5. 



V. CONCLUSION 



In conclusion we have studied the propagation of a bright matter wave soliton through 
a barrier potential with rapidly oscillating position. We have shown that the soliton can 
be dynamically trapped by such a potential. It would be interesting to extend this model 
to the case of oscillating barrier and standing well. In this case the effective potential will 
have the form of two barriers separated by a well. Applying the results obtained in Ref. [^] 
we can expect for the transmission of the nonlinear matter waves such phenomena as the 
instability and quantum turbulence on the time scales of the quantum dwell time in the 
well. Another physically relevant system is the rapidly oscillating barrier in the well. The 
effective potential in this case simulates a heterostructure and effective nonlinearities can 
exhibit chaotic dynamics in the tunneling process [lfll |. 

In the case of a well with rapidly oscillating position the effective potential will have a 
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double-well form. In such a configuration we can expect the existence of dynamically induced 
macroscopic quantum tunnelling and localization phenomena, i.e. dynamically generated 
bosonic Josephson junction. These problems require of a separate investigation. 
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